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Abstract 

We study a confined system of Dirac fermions in the presence of inhomogeneous magnetic 
field. Splitting the system into different regions, we determine their corresponding energy spectrum 
solutions. Wc underline their physical properties by considering the conservation energy where 
some interesting relations are obtained. These are used to discuss the reflexion and transmission 
coefficients for Dirac fermions and check the probability condition for different cases. We generalize 
the obtained results to a system with gap and make some analysis. After evaluating the current- 
carrying states, we analyze the Klein paradox and report interesting discussions. 
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1 Introduction 



The Dirac formalism plies an important role not only from mathematical point of view but physical 
one as well. The recent observation of the anomalous quantum Hall effect in graphene [H [2] renewed 
the interest to this formalism. In fact, many questions, raised in graphene, found their solutions by 
adopting the Dirac formalism as cornerstone. Among them, we cite the confinement [3] that is much 
needed to describe the transport properties in graphene. Subject that attracted much attention where 
interesting developments appeared by dealing with different issues, for instance we refer to [HE]. 

On the other hand, the quantum wires (electron waveguides) with quantized conductance can 
be formed in graphene [UJ. Such electron waveguides are indispensable parts of any conceivable all- 
graphene device. In lithographically formed graphene ribbons, the electronic bandstructure is theoret- 
ically expected to very sensitively depend on the width and on details of the boundary [7] . On top of 
that, disorder and structural inhomogeneity are substantial in real graphene [8]. For narrow graphene 
ribbons or electrostatically formed graphene wires [U], conventional conductance quantization thus 
seems unlikely [ID]. This expectation is in accordance with recent experiments 

Magnetic barrier technology is well developed [HI [13] Q3] and its application to graphene sam- 
ples appears to pose no fundamental problems [6j. In fact, snake states are experimentally studied 
in other materials [12\ I15j. mainly motivated by the quest for electrical rectification. On the theory 
side, the confined Schrodinger fermions in the magnetic field (with B' = 0) is discussed [16] as well 
as the asymmetric cases [IT]- For the Dirac- Weyl quasiparticles encountered in graphene, however, 
such calculations are not reported. The inhomogeneous magnetic field case in graphene is analyzed 
in |18j . Theoretically, the electron waveguides, in graphene created by suitable inhomogeneous mag- 
netic fields, is considered [6j. The properties of uni-directional snake states are discussed. For a certain 
magnetic field profile, two spatially separated counter-propagating snake states are formed, leading to 
conductance quantization insensitive to backscattering by impurities or irregularities of the magnetic 
field. 

Subsequently, The tunneling effect of two-dimensional Dirac fermions in a constant magnetic field 
is studied [19] , This is done by using the continuity equation at fixed points to determine the corre- 
sponding reflexion and transmission coefficients. For this, a system made of graphene, as superposition 
of two different regions where the second is characterized by an energy gap t', is considered. In fact, 
concrete systems are treated to practically give two illustrations: barrier and diode. For each case, 
the transmission in terms of the ratio of the energy conservation and t' are discussed. Moreover, the 
resonant tunneling by introducing a scalar Lorentz potential is analyzed where it is shown that a total 
transmission is possible. 

Motivated by the above progress and in particular [6j [19] , we deal with other features of the system 
considered in [6j. Such system is composed of different regions submitted to two magnetic fields and 
confined to a constant potential. This allows us to treat each region separately by determining the 
corresponding energy spectrum solutions. We underline some physical properties of their spectrum 
by taking into account of the energy conservation where interesting relations are obtained. Using 
the continuity at different points, we explicitly evaluate the reflexion and transmission coefficients. 
Combining all, we chow that the probability condition is well verified. As second task, we consider the 
present system with energy gap t' and do the same job to derive its eigenspinors as well as eigenvalues. 
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It is shown that even the reflexion and transmission coefficients take new forms in terms of gap but 
the probability condition still verified. Interesting limits are discussed, which concern total reflexion 
and transmission of the system with gap. 

Finally, we treat the Klein paradox by using the current-carrying states where different limits 
and discussions are presented. More precisely, we evaluate the currents for each region and use their 
relations to the reflexion and transmission coefficients to check the probabilities. Subsequently, three 
different cases are considered, which correspond to week, intermediate and strong potentials. We 
notice that two last cases are shown negative transmissions. However, by combining all coefficients we 
end up with a sum equal unity. 

The present paper is organized as follows. In section 2, we consider a confined Dirac fermion in 
inhomogeneous magnetic field ([T]). After getting the eigenvalues and eigenfunctions, we analyze the 
energy conservation that allows us to derive interesting relations between involved quantum numbers 
and parameters. In section 3, we study scattering between two regions to determine the reflexion and 
transmission coefficients, which will be used to discuss the probability conditions of the present system. 
We do the same job in section 4 but by considering three regions where the first is equivalent to the 
third. The continuity at each point leads to express the coefficients entering in the game in terms of 
different parameters. In section 5, we introduce a gap like a mass term and analyze the tunneling 
effect of such case. We study the Klein paradox in section 6 by involving the currents corresponding 
to different regions and consider three cases. Finally, we close by concluding our work. 



2 Dirac fermions in inhomogeneous magnetic field 

We consider a system of massless Dirac fermions through a strip of graphene characterized by the 
length d and width W in the presence of inhomogeneous magnetic field. More precisely, we introduce 
two magnetic fields B and B' , such as 





' B, 


x < —d 


B(x) = < 


B', 


\x\ < d 






x > d. 



According to the configuration ([T]), we decompose the present system into three regions. Schematically, 
we end up with Figure 1 
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Figure 1: Magnetic field profile. 



Clearly, regions I and III are similar but different with respect to region II. Note that, the system 
characterized by Figure 1 has been analyzed in [6] for possible quantum wires in graphene. However, 
in the present work we study other features of such system to deal with different issues, which concern 
tunneling effect and Klein paradox. 
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2.1 Dirac Hamiltonian 



Before writing down the appropriate Hamiltonian of the system (Figure 1), let us derive the corre- 
sponding gauge field to the configuration (pTJ) . Indeed, using the continuity of the potential to obtain 

A\(x) = Bx+(B- B')d, x < -d 



A } (x) 



\x\ < d 
x > d. 



(2) 



A\\{x) = B'x, 
A\u(x) = Bx-(B- B')d, 

where j is labeling regions 1,11 and III . It clear that, for B = B' we end up with one potential and 
therefore three regions become similar to each others. 

In the systems made of graphene, the two Fermi points, each with a two-fold band degeneracy, can 
be described by a low-energy continuum approximation with a four-component envelope wavefunction 
whose components are labeled by a Fermi-point pseudospin = ±1 and a sublattice forming an hon- 
eycomb. Specifically, the Hamiltonian for one-pseudospin component for the present system can be 
written as 



iTj = V F U ■ 7Tj + Vj(x) 

where the components of the conjugate momentum tt- } = p + -Aj are given by 

e 

Px, TTyj = Py + -Aj(X) 



7T 



(3) 



(4) 



and Vj(x) is the potential barrier that has a rectangular shape, which is infinite along the y-axis and 
has the form 



V] {x) = { "°' — (5) 

I 0, otherwise 

where Vq > 0. Injecting all in ([3]) to get 

Vj(x)/v F p x - ip y - if A- } (x) \ 
Px + ipy + i|Aj(a?) Vj(x)/v F J 

At this stage, it is convenient to introduce the annihilation and creation operators. They can be 
defined as 



(6) 



Oj = iPx+Py + -Aj(x), 

which obey the canonical commutation relations 

' 2h 2 



-iPx +Py + -Aj(x) 
c 



(7) 



aj,a: 



l \B'\ 

2h 2 



II 
III 



(8) 



where the magnetic lengths lg = y ^ and l\ B i\ - 

and B', respectively. The Hamiltonian (0) can be written in terms of aj and aj as 



^g7| are corresponding to the magnetic fields B 

t 



H\ = ivp 



V s (x)/iv F 



,t 



(9) 



Vj(x)/iv F J 

which is encoding all regions. This will be used to study each region separately and derive the 
corresponding energy spectrum solutions. 
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2.2 Energy spectrum solutions 

We determine the eigenvalues and eigenspinors of the Hamiltonian Hy Indeed, the the Dirac Hamil- 
tonian describing region I is obtained from @ as 

7). do 



The operators a\ and a\ can be rescaled to define others, such as 

V2h ah 1 ~ V2h 



5 i = -J^ a ^ b \ = -7^Z a \ ( U ) 



which verify 



Using these to write H\ as 



[b u b\)=l. (12) 



m = ihu c y ° o 6 ' j (13) 

where we have set uo c = V^t^ 1 as a cyclotron frequency. 

To get the energy spectrum solutions of (|13|> . we solve the eigenvalue equation for a given spinor 
I <Pi 



<P2 



of H\. This is 



which is equivalent to 



bA" 1 (14) 



To 1 : - : 



and leads to two relations between spinor components 

— iftw c b\ip2 = E\tp\ (16) 
ifajj c b\(pi = E{ip2- (17) 

Now inserting ()16|) into (|17p to obtain a differential equation of second order for if2 

h 2 u 2 c b}b m = Efip 2 . (18) 

It is clear that ip 2 is an eigenstate of the number operator n = b\b\ and therefore we identify ip2 to the 
eigenstates of the harmonic oscillator | n), namely 

tf2 H n) (19) 

and its eigenvalues read as 

E\ !n = shuj c ^/\n\ (20) 
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where n is obviously the eigenvalues of n, with n = 0, ±1, ±2, • • •, and s = sgn(re). Note that, n = 
corresponds to the lowest Landau level, i.e. zero mode energy. 
Using (fl~6|) . (fl~9j) and (f20|) to get the first component as 



E\ 



: k | n) 



-is n — 1) 



(21) 



which gives the eigenspinors 



0l,n 



-si | n — 1) 
I 



(22) 



In terms of the parabolic cylinder functions D n (x + xqi) [20], the eigenspinors in the plane (x,y) are 



<P\,n,k y (x,y) 



where xq\ = k y l 2 B + ( 1 — ) d and D n (x) are related to Hermite polynomials via 



1 / -^|„|_i(x + x i) 
V% I D| n | (x + ccoi) 



(23) 



D n (x) = (/ B 0Fn!2 n ) 2 exp 



fln(x). 



(24) 



As far as region II is concerned, we use the mapping n — > m and s — > s' in (j23|) to obtain the 



eigenspinors (f>\\,m,k y {x,y) in terms of X02 = Qy^fs'] an d * ne corresponding eigenvalues 

-Ell,; 



s hw r \/ \m\ + Vr 



(25) 



where s' = sgn(m) and uj' c = is the cyclotron frequency associated to the magnetic field B'. 

Finally, for region III the eigenvalues and the eigenspinors are similar to those of region I except 
that the correspondence 



xo3 = k y l 2 



y l B 



~B 



d 



(26) 



must be taken into account in (|23p . Note that, for B' = —B (with B' < 0) the eigenspinors for three 
regions can be expressed with the same position xq = k y lg. 

2.3 Illustrations 

To give some illustrations, we focus on the eigenfunctions of four lowest states in region II, which are 
summarized in Table 1. Different plots are given in Figure 2 those show that the m th eigenfunction 
has m-nodes, namely there are m-values of x for which 0n jm (x) = 0. 



Number 



Energy eigenvalue 



Energy eigenfunction 



m = 



ii,o — 



<£i i,o 0*0 



e -( x+X02 f/2l* Bfl 



m = 1 



E ni = hw' c + V 



0n,iO*O 



i 



2l\B' I yfx 



2 2 ( S+202 ) e -{x+X02?/2l 2 ]Bl] 



l \B'\ ) 



m = 2 



Eu 2 = yftHut + Vo 



A I X+XQ2 



-{x+x 02 ) 2 /2l 2 . 



m = 3 



E n ,3 = V3huj' c + V 



0ii,3(^) 



X+XQ2 



B< 



12 



X + XQ2 
l\B>\ 



e -(x+x 02 ) 2 /2l 2 B , 
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Table 1: Normalized eigenfunctions for four lowest states of a one-dimensional potential energy field. 






Figure 2: The spatial shapes of the eigenfunctions 4>\\ t m(x), with m = 0, 1, 2, 3 and l\g 
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On the other hand, the wavefunctions have observable properties. Indeed, if the position coordinate 
is changed from x to -x, the eigenfunction has a definite symmetry 

( 4>\\,m(-x) = +<j)\\, m (x), if m is even 

\ <t>\\,m{-x) = -(f>\\,m,(x), if m is odd 
which is nothing but the parity symmetry. Furthermore, the position probability density of fermion 
is given by 

\<l> m (x,y)\ 2 = \<l> m (x)\ 2 . (28) 

By plotting this for some specific values of m, we deduce an interesting conclusion. From Figure 3, 
we notice that the fermion can have any location between x = — oo and x = +oo, in marked contrast 
with a classical fermion, which is confined to the region —A < x < +A where A is the amplitude of 
oscillation. 



« 0.25 - 



J - 
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Figure 3: The position probability densities |</>ii )m (:c)| 2 corresponding to Table 1. 
2.4 Energy conservation 

In the interface between regions, there is conservation of the tangent components of the wave vector, 
i.e. k y = q y , and conservation of the energy. This is 

E\ = E\\ = E (29) 

which leads to the constraint 



\n\ \B'\ E 2 
\m\ B (E- V ) 



2- (30) 



Since n and m are integer values, the r.h.s term must be a fractional number, which can be written as 

\n\=K\m\, KeQ + . (31) 

This relation is very important because without such set one can not talk about tunneling effect in 
the present case. We will clarify this statement from next section and exactly when we begin by 
calculating different quantities in order to check the probability condition. 
Now let us return to (1201) and (1251) to write the ratio as 



* = (32) 



n\ — 1 -g- L |?Ti| 



Recall that in the region II we have Vq > 0, which implies that the energy E can be either positive or 
negative. Therefore, we should distinguish between two situations as listed below 



E > (s = s' = +1) 


E < (s = s' = -1) 


i i ^ \ B '\ i i 


i i / \ B '\\ i 
n < J -g- L to 



Table 2: Positive and negative energies and their corresponding quantum number configurations. 



These energies can be plotted to explicitly illustrate their behavior in terms of different quantities 
entering in the game, which are given in Figure 4. 
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Figure 4: Variation of the ratio p| in terms of ^-jj |?7i| for E > and r^r|^| for E < 



3 Two regions in inhomogeneous magnetic field 

To treat a concrete example of the present system, we consider a barrier submitted to an inhomoge- 
neous magnetic field. This barrier can be seen as superposition of two regions separated by an interface 
localized at a fixed point. We study the tunneling effect by evaluating the reflexion and transmission 
coefficients at interface, which in our case corresponds to the point zero. These will be used to show 
that the probability condition is exactly one and emphasis what makes difference with respect to 
without confinement case [19]. To do this task, we follow the same lines as has been done [19] and 
distinguish between propagation with positive and negative incidences. In both cases we deal with 
propagation from region I to II, II to III and vice verse. 



3.1 Propagation with positive incidence 

To proceed, let us first define the eigenspinors for three regions in positive and negative direction of 
the variable x. In region I, we write 



= J_ / -siD H _±(x + x 01 ) 
s/2 \ D\ n \(x + xqi) 

Similarly, in region II we have 

if -s'iD\ m \_ 1 (x + xo2) 
(p\\ + = —= 1 1 

V2 y D\ m \(x + x 2) 
as well as in region III 



<P\\ 



1 

71 



i 

71 



D\ n \(-X - Xqi) 



s'iD\m\-i(-% - X02) \ j k 

D\ m \(-X-X 02 ) j 



(33) 



(34) 



'Ml, 



1 / -siD^^x + X03) 
V 7 ^ y D\ n \(x + x 03 ) 



Jk y y 



5|||, 



1 
71 



-siD| nh i(-x-x 03 ) \ e i kyV 

D\n\{-X-X z) J 



(35) 



where (±) refer to the positive and negative propagations, respectively. 

We start by analyzing the case of propagation from region I to region II. Indeed, at the interface 
x = (for all y) , the continuity of the system gives 



^l,+ + r nm4>\ - = t n m4>\\, 



(36) 



s 



where r+ m and t^ m are reflexion and transmission coefficients, respectively, in positive propagation. 
From (1331) and AMI), we obtain 



/ — ai-D|„|_i(a;oi) ] + r + [ -s^|n|-i(-^oi) \ =t + ( ~ s> iD \m\-i{xQ 2 ) \ ^ 
\ D\ n \{x m ) J nm \ D\ n \{-x m ) J nm \ D| m |(cc 02 ) / 

They can be solved to get the coefficients in terms of some constants, which are magnetic field depen- 
dent. They are 

sA\ + s B\ 

t nm( x 01> x 02) = ~ j „ (39) 

where we have set 

A\ = A nm (x i,x 02 ) = D\ n \_ 1 (x 01 )D\ m \(xo2), 
B\ = B nm (xoi , x 2 ) = -D|m|- 1(^02 )D\ n \(xoi) 
C\ = C n (x i) = D| n |_i(x i)£>|„|(a;oi). 

On the other hand, considering propagation from region II to region I, the continuity at point zero 
reads as 

0N,+ + r tm$\\ - = ttin<t>\,+ (40) 

which implies 

/ s'iD\ m \_ 1 (x Q2 ) \ +r + ( s ' D \m\-i(-xoi) \ =t + I s£>| n |_i(x l) \ ^ 
\ D\ m \(x 2) J mn \ D\ m \(-x Q2 ) J mn \ D\ n \(x m ) J 

These lead to the solution 

r+»(*oi,*02) = (-l) |m| S 'f\~^\ (42) 



sA, + s'B ] + 



2s'F\ 
sA\ + s B\ 



where F\ is 



F\ = F m (x 02 ) = D\ m \_ 1 (xQ2)D\ m \(xQ 2 )- (44) 

In similar way, we show that the reflexion and transmission coefficients corresponding to propaga- 
tion from II to III are given by 

r+^sos) = (-1)'™' *?"" S 1W (45) 

sA\\ + s'B\\ 

C^02,x 03 ) = (46) 

SA|| + S D\\ 

where the constants read as 

■An = A nm (x 02 ,xo 3 ), B u = B nm (x 02 ,x 03 ). (47) 
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As far as the propagation from III to II is concerned, we find 

rnm{X02,X03) = (-1 1 1 48 

sA\\ + s'B\\ 

t+ m {x m ,xvz) = 2 ^'' (49) 

SA\\ + S D\\ 

where C\\ is 

C,i = C n (x 03 ). (50) 

This summarizes our analysis for propagation with positive incidence, which together will be used to 
discuss different issues and before doing so, we need to analyze negative incidence. 

3.2 Propagation with negative incidence 

We determine the reflection and transmission coefficients for the propagation with negative incidence, 
which will be denoted as r~ , and t~ , for the cases of the propagations from I to II, II to III and 

nm/mn nm/mn r r o i 

vice verse. To reply this inquiry, we use the same analysis as before but one should take into account 
the negative sign of variable. 

In doing our task, for region I we write the corresponding eigenspinors as 



1 / -siD {nhl (-x - x i) \ eikyy - 1 I - siD N-i( x + x oi) 

\/2 y D\ n \(-x-x i) J \/2\ D\ n \{x + x m ) 

For region II, we have 



(51) 



= 1 ( s^D lmhl (-x - x 02 ) \ ^ ^_ = J_( s'iD^x + xv) \ (52) 
V 2 V D \m\(-x-x 02 ) J V2 y D\ m \(x + X 02 ) J 

In region III, we write 

0m + = J_( siD^i-x - xos) \ ^ ^ _ 1 / siD^x + xos) \ eikyy _ m 
V 2 V D\ n \(-x-x 03 ) J V2 y D\ n \(x + x 03 ) J 

Now let us treat each case to evaluate reflexion and transmission coefficients. Considering the 
propagation: I — > II to obtain 

01,+ + r nm4>\- = tnmh,+ ( 54 ) 

at point zero. Injecting (f5"Tj) and ([52]) into (f54"|) . one gets 



r- m (x i,x 02 ) = (-l) 1 " 1 ^' (55) 

sA\ + s B\ 

t- m (x 01 ,xo 2 ) = (-1)IH+IH ^ 2 ^ , (56) 



For II — > I a similar relation to (|54|) reads as 

011,+ + r mn0N - = Km<h,+- ( 57 ) 
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The solutions are given by 

r- n (x 01 ,x 02 ) = (-l) H fjj'~*flj (58) 

^(wb) = (-l)IH + IH_j£^_. (59) 
After applying the same technique as before, we show that the propagation: 1 1 — > 1 1 1 gives 

r mn (xo2,x 03 ) = (-1)' l sAu + s , Bu (60) 

t^ n (x 02 ,x 03 ) = (-l)l"l+H 2S ' F \ (61) 

sA\\ + s'i?n 



and III — > II leads 



r nm (. 02 ,* 03 ) = (-1)1 + (62) 

(nm^oa.ato) = (-1)N+H /f" . (63) 

By inspecting the forms of different coefficients obtained so far, one can establish a symmetry 
between them. Indeed, We show the relation 

gjO^fog) = r|(soi.SQ2) = | B | + | TO | 
*| (^01,^02) rf i {x i,x 02 ) 

where the pair of index is chosen such as (i 7^ j) G {n, m}. Note that, the same relations are also valid 
for the couple (xo2i^03)- 

3.3 Reflexion and transmission amplitudes 

Now let us collect the products of our results by checking their importance. In fact, we discuss the 
reflexion and transmission amplitudes between regions to emphasis the influence of each parameter 
on them. For propagation between I and II, we define the reflexion and transmission amplitudes as 

p(x i,x 02 ) =rf j (x i,xo 2 )r^(x i,x 02 ), t(x 01 ,x 02 ) = tf j (x i,x 02 )tf i (xoi,Xo 2 ). (65) 
After replacing, we end up with 

[sAj - s'B{\ 2 

p{x i,x 02 ) = (66) 
[sA\ + s'B\\ 

T{x i,x 02 ) = — 7—^2=-—. 7— ^2 67 

[sA\ + s'B\} z [sA\ + s'B\Y 

where the relation C\F\ = A\B\ is satisfied. A straightforward calculation shows that the probability 
sums to unity, namely 

p(x i,x 02 ) + t(x i,x 02 ) = 1. (68) 

To stress how the amplitudes behave in terms of different parameters, we give Figure 5. It is clear 
that, p(xqi, X02) and t{xqi,xq 2 ) change with respect to variation the magnetic field B' for different 
values of B, m, d and k y . 
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Figure 5: Reflexion p(xoi,xo2) (red line) and transmission t{xq\,xq2) (green line) coefficients between 
region I and 1 1 for a magnetic barrier of width 2d at various magnetic field B' for two cases: a) B = 10, 
b) B = 15 and for different values of (m, d, k y ). 



It worthwhile to invert the situation by varying B for two values of B' a given configuration of the 
parameters (m,d,k y ). This summarized as follows 
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Figure 6: Reflexion p(xqi,xq2) (red line) and transmission t(xoi,^o2) (green line) coefficients between 
region I and II for a magnetic barrier of width 2d at various magnetic field B for two cases: a) B' = 10, 
b) B' = 15 and for different values of (m, d, k y ). 



As far as the propagation between II and III is concerned, we use the same definition as above to 
write the amplitudes 

P(x 02 ,X 03 ) =r^(x02,S03)7^(X02,Z03), t(xq 2 ,X Q?) ) =tf j (xQ2,XQz)tf i (x 2,Xo ?l ) (69) 

which give 

p{x 02 ,xw) = —J—T2 ( 7 °) 

[sA\\ + s'B u \ 

Ass'C\\F\ iss'AuBu 

T{X02, X 03 ) = = 7T (71) 

[sAn + s'Bn] 2 [sAn + s'Bu] 2 

where C\\F\ = A\\B\\. Using these to verify 

P{XQ2 ,X 03 ) + T (x 2 , X 3 ) = 1 • (72) 

As before we illustrate this result by making different plots, which are given by Figure 7 
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Figure 7: Reflexion p(xq2,xq^) (red line) and transmission t(xo2,^03) (green line) coefficients between 
region II and III for a magnetic barrier of width 2d at various magnetic field B' for two cases: a) 
B = 10, b) B = 15 and for different values of (m, d, k y ). 
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Figure 8: Reflexion p{xq 2 ,x 03 ) (red line) and transmission t(xq 2 ,xq 3 ) (green line) coefficients between 
region II and III for a magnetic barrier of width 2d at various magnetic field B for two cases: a) 
B' = 10, b) B' = 15 and for different values of (m, d, k y ). 

These are among the interesting results derived so far. In fact, it tell us the transmission of barrier 
in inhomogeneous magnetic fields can not be greater than one, which is analogue to what obtained in 
one field case and without confinement |19j . 

Let us present some discussions and derive interesting results those have applications in physics 
areas. Indeed, similar relations to what obtained above exist also for the photon optics cases. For 
instance, one can write ([68]) and ([72]) when a light beam is reflected and refracted in a diopter between 
regions I- 1 1 and 11-11 1. The first case is characterized by the configuration 



For n = 



p{xqx,xq 2 ) = 1, t(x 01 ,x 02 ) = 

P(xq 2 , Xqz) = 1, t(x 02 , X 03 ) = 0. 



This is an expected result since the transmission r = 1 — p must be zero in the case where the wave in 
the n-region enters in the m-region and vice verse. In fact, the interface between tree regions behaves 
like a mirror where the reflexion is total. The second case is described by 



For s = s' and n = ±m 



p(xoi,x 02 ) = 0, t(x i,x 02 ) = 1 
p{xo2, x 03 ) = 0, t(x 02 , x 03 ) = 1 



which means that the interface between region l-ll and 11-11 1 behaves like a non-reflective diopter, 
namely there is a total transmission. Note that, these two cases have interesting interpretation in 
optics physics. 



4 Tree regions in inhomogeneous magnetic fields 

We study another case of a physical system composed of a region indexed by the quantum number m 
of length 2w separating two others indexed by the same quantum number n. This will allow us to see 
how the above results will be changed to the present case and underline what make difference with 
respect to the former analysis. 

4.1 Reflexion and transmission coefficients 

The present situation is quiet different from the former one. We use the above tool to write the 
continuity equation at the points x = —d and x = d. Then, we derive the quantities needed to discuss 
the reflexion and transmission coefficients as well as the corresponding probabilities. This will be done 
by treating propagation with positive and negative incidences. 

We study the positive incidence by considering the geometry that corresponds to the first interface, 
which is 



15 



D 





i 


-d 


B' 


X 



Figure 9: Magnetic field profile between regions I and II. 

Using the above spinors to setup the continuity equation at the point x = —d. Doing this process to 
get the relation 

((>\,+(-d) + = a^u, + (- d ) + PHA-d) ( 73 ) 

where r + is the reflection coefficient for positive incidence, which will be determined together with the 
parameters a and /3. ([73]) gives 

(74) 
(75) 



S-D| n |-i(^l) + r + sD\ n \A-d\) = as' D\ m \_ t {d 2 ) + fts' D\ m \_ x (-d 2 ) 
D\ n \{di) + r + D\ n \(-d{) = aD\ m \(d 2 ) + PD\ m \(-d 2 ) 

with the constants d\ = xoi — d and d 2 = £02 — d. These can be solved for a and /3 to obtain 
sA nm (d 1 ,d 2 ) + s'B nm (d 1 ,d 2 ) + r+(-l)l n l [s'B nm (d 1 ,d 2 ) - sA nm (d 1 ,d 2 )] 



a 







2s>F m (d 2 ) 

s'B nm (di,d 2 ) - sA nm (di,d 2 ) + r+(-l)l n l [sA nm (di,d 2 ) + s'B nm (di,d 2 )] 



2(-l)\m\ s 'F m (d 2 ) 
In terms of the reflexion and transmission coefficients, we have 



1 



O: 



tmn(dl, d 2 ) 
(_l)l n l + l m l 



- r nA d i,d 2 ) 

tmn(di,d 2 ) 



r nm(dl,d 2 ) 



tmn(di,d 2 ) t mn (di,d 2 ) 
To accomplish such analysis we consider the second interface as shown below 

B 



(76) 
(77) 

(78) 
(79) 



1 


d 




B' 




► 

X 


Figure 10: Magnetic field profile 


between reg 


we have 






a4>\\A d ) +P<t>\ 


Ad) 


= t + ^\\A d ) 



(80) 

where t + is transmission coefficient for positive incidence. After replacing, we end up with a system 
of equations, such as 



as'-D| m |_i(4) +/3s / £>| m |_i(-4) = * +s ^|n|-i(^3) 
aD H (d' 2 ) +/3£> H (-4) = t+D ]n] (d 3 ) 



(81) 
(82) 
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with d' 2 = d + xq2 and ^3 



a 



d + X03. The solution reads as 

ran 

{d 3 ,d' 2 ) + s'B nm (d 3 ,d' 2 ) 



2s'F m (d' 2 ) 



r 



tmn(d' 2 , d 3 ) 



(83) 



p = t+(-l)M 



s'B nm (d 3 ,d 2 ) - sA nm (d 3 ,d' 2 ) 



h+h: 



tmn{d 2 , d 3 ) 



To determine the coefficients for positive incidence we simply use (| 7811 79 j) and (|83H84p . Combing 
obtain 

„+ r nm (di , d 2 ) - r~ m (d' 2 , d 3 ) 



1 - r nm (di, d 2 )r nm (d 2 , d 3 ) 



(84) 
all to 

(85) 



tmn(di, d 2 ) 
tmn(d\, d 2 ) 



1 ~ r nm{di,d 2 )r nm (di,d 2 ) 
1 - rnm(di,d 2 )rnm(d 2 , d 3 ) 
T(di,d 2 ) 

1 



(di,c? 2 )r nm (ii / 2 ,d3) 

Now let see how the above results will be written by considering the negative incidence 
Indeed, applying the same machinery as before to get 

r Xm (di , d 2 ) - r + (d' 2 , d 3 ) 

v = — 1 ■ 

1 - rt,m{d\, d 2 )rt,m(d 2 , d 3 ) 

r n,m(di,d 2 ) - r- m (d' 2 ,d 3 ) 
1 - rn,m(d\, d 2 )rn,m(d 2 , d 3 ) 



case. 



(87) 



tn,m(d' 2 , d 3 ) 
tn.m(dl,d 2 ) 
tn,m(d' 2 , d 3 ) 
tn,m(di,d 2 ) 



r(di,d 2 



1 - rZ, m (di,d 2 )rn,m (d 2 ,d 3 ) 

r(d 1 ,d 2 ) 

1 



(di,d 2 )r n , m (d' 2 ,d 3 

Having obtained the above results, we analyze the corresponding probability and give comments. This 
issue and related matter will be considered in the forthcoming subsection. 

4.2 Probability 

To characterize the behavior of the present system, we study the incident beam. This can be achieved 
by calculating the probability of reflexing and transmitting beam. Indeed, let us adopt the definition 



R 



r r 



T = t + t~ 



After calculation, we find 



R 



p(di,d 2 ) + p{d' 2 ,d 3 ) - 2r- m (di,d 2 )r- m {d' 2 , d 3 ) 
1 + p(di,d 2 )p(d 2 ,d 3 ) - 2rn, m (di, d 2 )rn,m(d 2 , d 3 ) 

1 + p(di,d 2 )p(d' 2 ,d 3 ) - p{di,d 2 ) - p(d' 2 ,d 3 ) 
1 + p(di,d 2 )p(d' 2 ,d 3 ) - 2r~,m.(di,d 2 )r~ :m (d' 2 ,d 3 ) 



(89) 

(90) 
(91) 
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Combining all to end up with probability 

R + T = l. (92) 
From this, we summarize the following conclusions: 

• The probabilities of reflection and transmission sum to unity, as must be the case, since they 
are the only possible outcomes for a fermion incident on the barrier. 



• (|90j) and (|91j) yield that under resonance conditions: 

|n| = \m\, s = s' 
the barrier becomes transparent, i.e. T = 1. 

• More significantly, however, the barrier remains always perfectly transparent for \n\ = \m\. 

• T = 1 is the feature unique to massless Dirac fermions. 

• T = 1 is directly related to the Klein paradox in quantum electrodynamics. 

5 Introducing gap 

In the present study, we consider the confined system in inhomogeneous magnetic field given by the 
configuration (1) but in the presence of an energy gap t' in the region II. We will show how the above 
results will be generalized to the gap case. 

5.1 Hamiltonian formalism 

As far as regions I and 1 1 1 are concerned, the eigenvalues and the eigenfunctions are those given before 
for case of without t' . However, in region II the Dirac Hamiltonian can be written as 

m i = < = v F aTT + V + t'a z . (93) 

Clearly, the mass term t'a z makes difference with respect to the former analysis. In fact, it will play 
a crucial role and lead to discover interesting results. In terms of matrix, H\\ takes the form 

/ Wvfc + tf \ 

V Vx + iPy + KA 2 (x) I \ V„ - f I 



For next purpose, we determine the energy spectrum solutions of (|94p . In doing so, let us fix 
(>n = I * J as a spinor of H\\ in presence of an energy gap t' , such as 



V n 2 



which implies two relations 



-ihuj' c D 2 v' 2 = (E u -Vq- t')^ (96) 
ihu'cD+rt = (E U -V + t)(f/ 2 . (97) 
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They are showing 

h 2 oj' c 2 D+D 2 <p' 2 = [{E„ + V ) 2 - t' 2 ] <p' 2 . 
It solution gives the second spinor component as 

V>2 = D \m\ (x + x 02 ) e iqyV , me; 
From (|98p . it is easy to obtain the eigenvalues 

E n , m = s'^h 2 oj' c 2 \m\+t' 2 + V Q . 



r =o, i'=o 

\m\ - 1 \m\ -2 M = 3 | 



Energy E 




(98) 



(99) 



(100) 



Figure 11: Landau levels E m in terms of magnetic field B' for different values of m. 
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Figure 12: Energy EV^ for a magnetic barrier of width 2d: a) as function of B' and t', b) E m as 
function of i?' and Vq. 

We notice that in the presence of energy gap, there is a gap separating the conduction and valence 
bands, which is missing in the case where t' = (Figure 12b). 

To complete the derivation of eigenspinors, we determine the first component ip^. Then, from (|96p 
and (1971) we obtain 



E\\ -V -t 



7 D| m |_i(x + a;o2)e 



iq y y 



(101) 
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After normalization the eigenspinors read as 

H mq (x,y) = ± ( -«^IH-i(- + -o 2 ) | w (1Q2) 
V2\ b m D lm] (x + x 02 ) 



where the constants are given by 



, En-Vp + M E U -V - s't' 
a m = s W — , b m = W — . (103) 

As concerning regions I and III, the corresponding eigenspinors 4>\, n ,k y {x, y, xoi) and 4>\\\,n,k y {x-, V, ^03) 
can be written in compact form as <fi\ t n,k y (x,y,xo) where 

x = kyl% - (l - x (104) 

such that x = —d and x = d give xo = x$\ and x = X03, respectively. 

On the other hand, the energy conservation between regions I and II gives 

E\ = E\\ = E. (105) 



After replacing, we show that the quantum numbers n and m verify the relation 

\n\ _ \B'\ E 2 

H ~ B (E-Vof-t' 2 ' 



(106) 



We have some remarks in order. In region 1=111 we have Vq = and tf = 0, thus (|106p reduces to 

IB'I 



\m\ B 

However, in region II there two cases: 



(E - Vq) 2 > t' 2 =^ \m\ = +m 
(E - Vq) 2 < t' 2 |m| = -m. 



Finally, the analogue of (132j) is now given by 

E J\n\ 



(107) 



(108) 



(109) 



5.2 Reflexion and transmission coefficients in the presence of if 

We will see how the results obtained before can be generalized to the present case. To proceed, 
we consider two (barrier) and three regions (diode). For barrier, we show that the reflexion and 
transmission coefficients are 

\ [sb m A\ - a m B\} 2 

P(X0,X 02 ) = —2- (HO) 

[sb m A\ + a m B\\ 



t'(xo,xq 2 ) 



Asb m a m C\F\ 
[sb m A\ + a m B\] 2 
4:sb m a m A\B\ 



[sb m A\ + a m B\Y 
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Clearly, what makes difference with respect to the former results is the appearance of the constant 
parameters a m and b m . These coefficients can be used to verify the probability condition 

p'(x , x 02 ) + t / (x ,x 02 ) = 1. (112) 



a) 
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Figure 13: Reflexion p'(xo,xo2) (red line) and transmission t'(xq,xo2) (green line) coefficients for a 
magnetic barrier of width 2d at various energy E for two cases: a) Vo = 15, t' = 5 and b) Vq = 30, 
= 10. 



At (Vo = 15, t' = 5) and (Vo = 30, t' = 10) the transmission and reflexion profile (Figure 13) between 
regions I — II and II — III show 

• —t' < E <t'\ t'(xq,xq2) — > 1 and p'(xo,xo2) — > 

• t' < E <\q — t': t'(xq,xq2) decreases and p'(xo, X02) increases. 

• Vq — t' < E < Vo + 1': there is no transmission and non reflexion (not allowed states). 

• E > Vq + 1'\ t'(x , xq 2 ) 1 and p'(x , X02) -> 0. 

As far as three regions are concerned, one can inspire from the case without gap to obtain the 
reflection and transmission amplitudes, such as 

R , _ ftjdudy) + p'(d 3 , d' 2 ) - 2r+ m (d 1 ,d 2 )r+ m (d 3 , d' 2 ) 
1 + p'(di,d2)p'(d 3 ,d' 2 ) - 2r^ m (d 1 ,d 2 )r^ m (d 3 ,d 2 ) 



T , = l + p'(d 1 ,d 2 )p'(d' 2 ,d 3 ) - p'{d 1 ,d 2 ) - p'(d' 2 ,d 3 ) 
1 + p'(di,d 2 )p'(d 2 ,ds) - 2rt,m(di, d 2 )rt, m (d 2 , d 3 ) 



After a straightforward calculation, we find 

R'+ T'=l. (115) 



Note that, for B = B 1 we discover the results obtained in [19], which shows that our results are more 
generals. 
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6 Klein Paradox 



We complete the present work by analyzing the Klein paradox for the present system. This can 
be done by introduce other considerations based on the current-carrying states and study different 
limiting cases. 

6.1 Propagation from region I to region II: (x = —d) 

We consider the scattering of a Dirac fermion of energy E from an electrostatic step-function potential 
as shown in Figure 14. This problem is an archetype problem in nonrelativistic quantum mechanics. 
For relativistic quantum mechanics, we will find that the solution leads to a paradox (Klein paradox) 
when the potential is strong. 

V(x) 











to 

r 


- E 




region I 






region II 



r «1 



Figure 14: Electrostatic potential idealized with a sharp boundary, with an incident free scalar fermion 
of energy moving to the right in region I. 



According to the previous analysis for two regions, it is easy to note that the solution of the Dirac 
equation in the region I and II are given by 




-siD\ n \_ x {-x - x i) 
D\ n \(-x - xqi) 



(116) 
(117) 



1 / -si-D| n |_i(cc + X l) 
y/2 \ D\ n \{x + xoi) 

J_ T ( -a m iD\ m \_ 1 (x + x 2 
V% \ b m D\ m \(x + x 02 ) 

where R and T are reflected and transmitted coefficients, respectively. Imposing the boundary condi 
tion that 4> be continuous at (x = —d) gives the relation 

sD\n\^l(-d + X 01 ) \ +R f sD \n\~l(d - X i) 

D\ n \{d-x m ) 



D\ n \(-d + x 01 ) 



m .D| m |_i(-(i + x 02 ) 
D\ m \(-d + x 02 ) 



(118) 



Solving for R and T to obtain 



R 
T 



sb m uiv 3 - a m viU3 
a m v 2 u 3 - sb m u 2 v 3 

s(uiv 2 - V\U 2 ) 
a m v 2 u 3 - sb m u 2 v 3 



(119) 
(120) 



where we use the notation 



u x = D H _ l (-d + x 01 ), 

v l = D \n\{-d + X 01 ), 



u 2 = D\ n \-i{d-x 01 ), 

V2 = D\n\(.d-X 01 ), V 3 



u 3 = D\ m \-i(-d + x <n) 
: D\ m \(-d + x 02 ). 



(121) 
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To proceed further, we introduce the current-carrying states. This is based on the current associ- 
ated to Dirac equation, which is 

J = evJ2(j) + ai<f) (122) 

i 

where i = x,y. As an immediate application, the incident current is given by 

Ji = ev (4>^~<Jx<P\ + 4 ) t (T y ( t ) \) = sevu\V\. (123) 

This can be used to evaluate the final currents to the left and right of the potential boundary, which 
read as 

Jx<-d = sev (ui + Ru 2 ) (v\ + Rv 2 ) (124) 
J x >-d = ^-\T\ 2 u 3 v 3 [i(a m b m - a m b m ) + (a m b m + a m b m )] . (125) 

(I125D can be split into three parts 



Jx>- 



-ev\T\ 2 a m b m u 3 v 3 , Vq > E + t' 

-iev\T\ 2 a m b m u 3 v 3 , E - t' < V < E + 1' (126) 

k ev\T\ 2 a m b m u 3 v 3 , V < E - t' 



with the condition a m bm > 0. 

Recall that, the reflexion and transmission amplitudes are related to the currents through 

R= J| ~ Jx< -\ T = ^>-d_ (127) 
Ji J| 

After replacing, we end up with 

R _ UiVi {a m V2U3 - sb m u 2 v 3 ) 2 - amb m u 3 v 3 (uiv 2 - viu 2 ) 2 
u x vi (a m v 2 u 3 - sb m u 2 v 3 ) 2 



U 3 V 3 [ijambm ~ flm^) + (ambm + ambm)] (U1V2 ~ V\ U 2 ) 

2mvi {a m v 2 u 3 - sb m u 2 v 3 ) 2 



2 



(129) 



These show that the probability is 

R + T = l. (130) 

We can inspect (|126p further to derive other results. This can be achieved by considering three 
interesting cases. 

6.2 Limiting cases 

In region II there are three distinct cases, depending on the strength of the potential. This is shown 
by Figure 15: 
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Intermediate Potential Region 
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Intermediate Potential Region 

, E-f 

Weak Potential Region 



Figure 15: Energy level diagram for a fermion in region II. 

Let us analyze each case separately and underline its physical properties. Indeed, in the weak 
potential that corresponds to E—Vq > t' , we have a restriction on the quantum numbers and parameter 
constants, such as 

\m\ = +m, a\ n = a m , $ m = b m . 
In such case, the reflexion and transmission amplitudes are given by 

R _ UlVl {a m V 2 U3 ~ sb m U2V3) 2 ~ a m 6 m M3»3 {uiV 2 - VlU 2 ) 2 

uivi {a m v 2 u 3 - sb rn u 2 v 3 ) 2 
T _ a m b m U3V3 (uiv 2 - viu 2 ) 2 ^ 132 ^ 
uivi (a m v 2 U3 - sb m u 2 V3) 2 

They verify the probability condition ()130p . Thus the incident beam is partly reflected and partly 
transmitted. This is similar to the result obtained in nonrelativistic quantum mechanics. The last 
expression shows that the total probability is conserved. 

As far as the intermediate potential is concerned, i.e. \E — Vq\ < t', different quantities reduce as 

|l7i| — 77t, — Q"im — ^m- 

The corresponding amplitudes are 

R _ UlVl («^2 M 3 ~ h m U l V 3) 2 + iambmU3V3 ~ VlU 2 f 

u x v x (a 2 m v\u\ - b 2 m ulvl) 2 

T _ ia m bmU 3 V3 (uiv 2 - v\u 2 ) 2 
uivi {a^vlul + b 2 m ulvl) 

where the probabilities sum to unity, as must be the case, since reflection and transmission are the 
only possible outcomes for a fermion incident on the barrier. 
In the strong potential case, i.e. \E — Vq\ > t', we have 

\m\ = -m, aL = -a m , bL = —b m 



which is showing 



2 



uivi (a m b m U3V3 (uiv 2 - viu 2 ) + a m v 2 U3 - sb m u 2 v 3 

R = ^ g — > 1 ( 135 ) 

uivi (a m v 2 U3 + b m u 2 V3) 

a m b m u 3 V3 (uiv 2 - viu 2 f 

1 = g < (loo; 

uivi {a m v 2 U3 + b m u 2 V3) 
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The probability is still conserved, but only at the cost of a negative transmission amplitude and a 
reflection amplitude which exceeds unity. The strong potential appears to give rise to a paradox. 
There is no paradox if we consider that in the strong potential case the potential is strong enough to 
create particle-antiparticle pairs. The antiparticles are attracted by the potential and create a negative 
charged current moving to the right. This is the origin of the negative transmission amplitude, i.e. 
(I136p . The particles, on the other hand, are reflected from the barrier and combined with the incident 
particle beam (which is completely reflected) leading to a positively charged current, moving to the 
left and with magnitude greater than that of the incident beam. 

7 Conclusion 

We considered a system composed of different regions of Dirac fermions in the presence of an inhomoge- 
neous magnetic field and confining potential V(x) in one-direction. The energy spectrum solutions are 
obtained in terms of different parameters and quantum numbers for each regions. To underline some 
physical properties of the obtained solutions, we analyzed the energy conservation. This allowed us 
to establish interesting relations and therefore solve some issues related to reflexion and transmission 
of the system. 

More precisely, by considering our system as as barrier, we derived interesting results. In fact, 
using the continuity equation at different points we explicitly determined the reflexion and transmis- 
sion coefficients. These are used to define the corresponding amplitudes and therefore to show that 
their probabilities sum to unity. Different cases are treated, which concerned total reflecting and 
transmitting beams where they are interpreted as mirror or diopter systems. 

Subsequently, we focussed on three regions of two fixed points d and —d. Writing the continuity 
at each point, we derived different quantities those are needed to characterize the beam of the present 
system. Indeed, we reached the conclusions that the probabilities of reflection and transmission sum 
to unity, as must be the case, since they are the only possible outcomes for a fermion incident on the 
barrier. Furthermore, (|90p and (|9ip yielded that under resonance conditions: |n| = \m\ and s = s' , 
the barrier becomes transparent, i.e. T = 1. More significantly, however, the barrier remains always 
perfectly transparent for \n\ = \m\. This latter is the feature unique to massless Dirac fermions and 
directly related to the Klein paradox in quantum electrodynamics. 

Another interesting cases is analyzed, which concerned introducing a gap. After getting the energy 
spectrum solutions, we discussed different issues and among them the energy conservation. This 
allowed us to generalize the former analysis to gap case. As interesting results, we showed that the 
probabilities of reflecting and transmitting amplitudes sun to unity as well. Requiring that B = B' , 
we recovered the result obtained in |19| . 

Finally, we discussed the Klein paradox by involving the current-carrying states for different re- 
gions. Using their relations to the reflexion and transmission amplitudes, we checked the probability 
by evaluating different quantities. Moreover, we treated three different limiting cases, which concern 
week, intermediate and strong potentials. For two last cases, the transmission amplitude is obtained 
with a negative sign, however when it is added to the reflexion coefficient gives a sum equal unit. 
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